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ABSTRACT 

We consider the implications of recent infrared and radio observations of the nucleus of 
NGC 4258. There is no direct evidence that the nucleus has been steadily accreting on the 
viscous timescale of the outer rnasing disk, which is ^ 10 9 yr. Thus the mass accretion rate in 
the outer disk need not be the same as in the inner accretion flow where most of the gravitational 
binding energy is released. We show that an advection-dominated flow model with a transition 
radius of ~ (10 — 100)GM/c 2 (where M is the mass of the hole) and M ~ 10 -2 is 

consistent with the observed spectrum from radio to X-rays. We also show that a thin (flat or 
warped) disk can fit the observed fluxes outside the X-ray band. The X-rays can be explained 
by means of a corona in such a model, but the absence of radio emission from the location of 
the putative central black hole provides a serious constraint on the properties of the corona. A 
wide range of accretion rates, 10 -4 <; M <; ICR 2 M 0 yr -1 , can be made to fit the data, but the 
most “natural” models have ICG 3 <) M lO _2 M 0 yr _1 . Physical conditions in the observed 
VLBI jet features can also be related to conditions in the inner accretion flow. We conclude 
with a list of future observations that might help to constrain the accretion rate. 

Subject headings: accretion, accretion disks — black hole physics - radio jets 


1. INTRODUCTION 

Most measurements of the physical properties of 

the maser spots allows one to determine the central 
mass (Perrnstein (1997)). By measuring maser spot 
accelerations (Haschick & Baan (1990), Haschick et 


al. (1994), Greenhill et al. (1995c), Greenhill et al. 


(1995a), Herrnstein (1997)|), or proper motions (Her- 

portunitios for advancing our knowledge, and bo it is — 

rnstein (1997)) one can also directly determine a ge- 


worth investing some effort to understand these ob¬ 
jects in detail. The orbiting water maser spots in 
NGC 4258 provide one of these rare instances. 

Water maser emission was first detected in the nu¬ 
cleus of NGC 4258 by |Claussen et al. (1984) . It 
was soon realized that the masers might lie in a disk 
( |Claussen &; Lo (1986)|) . Nakai et al. (1993)| dis¬ 


covered high velocity maser emission, which |Wat- 
|son &; Wallin (1994)| quickly interpreted as a natu¬ 
ral outcome of maser emission from a disk. VLBI 
observations (|Greenhill et al. (1995b)|, |Miyoshi et 
(1995)| 


al. 


Moran et al. (1995), Herrnstein, Green 


hill, &; Moran (1996H) demonstrated that the maser 


spots are positioned on the sky as one would ex¬ 
pect for a nearly edge-on thin disk.—A model fit to 


The measurement of an accurate mass and dis¬ 
tance eliminate two of the principal uncertainties that 
plague efforts to understand the physics of active 
galactic nuclei (AGN). The mass of the central ob¬ 
ject in NGC 4258 (which we will freely refer to as a 
black hole) is (4±0.25) x 10' M 0 (Herrnstein (1997); 
ler). The distance is measured to be 7.3 ± 0.3 Mpc, 
with independent and consistent results from mea¬ 
surements of m aser spot acceleration and maser spot 
proper motion ( Herrnstein et al. (1997c) ). 

One of the other main uncertainties that plagues 
AGN studies is the accretion rate M. The accretion 
rate in NGC 4258 has been controversial, and there 


are a number of estimates in the literature: Neufeld & 
Maloney (1995)(hereafter NM), 7x 10 -6 a_i M 0 yr -1 . 


1 Also Isaac Newton Institute, 20 Clarkson Rd., Cambridge CB3 0EH, UK 
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where a = 0.1a_i ^ 1 is the usual dimension¬ 
less angular momentum diffusion coefficient of accre- 


the jet radial velocity. More detailed models of the 
jet suggest an age of at least 10 6 to 10 7 yr (e.g. can 


(198l )|); |Lasota et al. (1996) , 0.014a_iM Q yr 


“1 


tion disk theory QShakura fc Sunyaev (1973) , Pringlcj Albada fe van der Hulst (1982) , Martin et al. (1989) ). 


|Cao &: Jiang (1997) , 8 x 10 3 q_i Mq yr 1 ; [Maoz 
Sz M jKee (1997)1 7 x 10 -3 M Q yr -1 ; [Kumar (1997) 


1.5 x 10 -3 Mq yr. These estimates differ because of 
disagreements about the nature and efficiency of an¬ 
gular momentum transport in disks and, to a lesser 
degree, because of disagreements about the heating 
and cooling physics in various parts of the accretion 
flow. The purpose of this paper is to clarify some of 
these differences, propose observational tests that re¬ 
duce the number of acceptable models, and thereby 
narrow the range of uncertainty in M. 

We begin in §2 by asking what the masing disk 
might tell us about the accretion rate. In §3 we give 
a critical review of the pertinent data; in §4 we con¬ 
sider models for the inner accretion flow; in §5 we 
consider implications of observations of the jet. A 
discussion and summary of interesting observational 
and theoretical problems follows in §6. 

2. MASING DISK 

Most estimates of the accretion rate in NGC 4258 
are based on observations of the masing disk and 
the assumption that the disk is in a steady state 
( dM(r)/dr = 0) from the masing disk in to the black 
hole, where most of the gravitational binding energy 
is released ( NM, |Cao &; Jiang (1997) , Maoz &; Mc¬ 
Kee (1997)1, Kumar (1997)|) . But is the masing disk 


Thus the jet provides a record over only a small frac¬ 
tion of the masing disk viscous timescale. The jet is 
roughly continuous over its length, however (there are 
no lengthy gaps in the jets: see for example Figure 
2 of jCecil et al. (1995a) ), which suggests approxi¬ 
mately steady accretion over a time tj. For standard 
a disk models for the inner accretion flow we find 
t v = t j ~ 10 6 yr at r ~ 10 3 GM/c 2 (the precise value 
depends on the accretion rate, a, etc.). Thus only a 
small region of the accretion flow in NGC 4258 needs 
to be in a steady state. 

3. CRITICAL REVIEW OF DATA 
3.1. Nuclear Spectral Energy Distribution 


in a steady state? In the absence of significant exter 
nal torques (e.g. a magnetohydrodynamic wind) one 
expects a steady state to be achieved on the viscous 
timescale t v = (all) -1 x ( r/H ) 2 , if mass is steadily 
supplied to the disk from r ^ 0.26 pc. 

In the masing disk, t v = 2.6 x 

10 9 a_i” 1 (200K/T)(r/0.2pc) 1 / 2 yr. 2 This is a sig¬ 
nificant fraction of a Hubble time. We know of no 
direct evidence that mass has been steadily supplied 
to the nucleus of NGC 4258 , or to any active nucleus, 
on this timescale. Thus there is no direct support for 
the assumption that the disk is in a steady state. 

NGC 4258 contains an unusual, kpc-scale jet de¬ 
tected in optical li nes (Pourtes fc Cruvellier (1961) ), 
radio continuum ( van der Kruit et al. (1972) ), and 
X-rays QCecil et al. (1995a) ). If the jet is pro¬ 
duced in the immediate vicinity of the black hole, 
and the jet mass loss rate is closely coupled to the 
accretion rate, as is commonly believed, then it pro¬ 
vides a fossil record of past accretion. The char¬ 
acteristic timescale for the jet is tj = R/Vj = 
2.4 x 10 6 (.R/5kpc)(2000kmsec -1 /Vj) yr where Vj is 

2 An optically thin coronal component with T ~ 8000 K, as 
timescale. But in the NM model, at least, most of the mass flux 
of the masing disk. 


Chary & Hcckhn (1997) have recently reported J,H, 
and K band (1.25,1.65, 2.21^m) observations of NGC 
4258 in O.b^ seeing. Using the J band image as a base¬ 
line, they subtract a scaled J band image from the H 
and K images and find a nuclear emission excess con¬ 
sistent with a point source. The excess is 1.1 mJy 
at H and 4.5 mJy at K. It is then possible to deduce 
the extinction and intrinsic luminosity of the source if 
the extinction follows a standard form deduced from 
the local interstellar medium ( Cardelli et al. (1989) ) 
and if the intrinsic spectrum has an approximately 
power-law form over the region of interest (/„ ~ v s )\ 

A v = 20.1(1 + 0.21s)mag, (1) 

log (vL v {H)/ ergs -1 ) = 41.63 + 0.32s, (2) 

\og{v L U {K)/ ergs -1 ) = 41.51 + 0.19s. (3) 

The spectral index s might reasonably range from 
~ 1/3 for the classical thin disk spectrum to ~ —1, 
which is a typical optical spectral index for an unob¬ 
scured AGN. The implied flux at J is 0.2(2.62) s mJy, 
which is small, justifying the use of the J band image 
as a baseline. Taking Ry = 3.1 ( [Spitzer (1978)1) , the 
extinction implies an obscuring column 

N h = 3.8 x 10 22 (1 + O.21s)(Z 0 /Z) cm -2 (4) 

(Z = metallicity) or, assuming cosmic abundance of 
helium, 


X = 0.09(1 + 0.21s)(Zq/Z) gem -2 . 


(5) 


These estimates are sensitive to the assumed ex¬ 
tinction curve, since 

dl„,L (fl) 53 _ 11S| (6) 

d\n(A H / Ak) 

suggested by NM, might carry a large M with a small viscous 
is carried in the cool, optically thick component at the inner edge 
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and 


dAy 


dln(Aff/X k ) 


= —50 — 1 D.s. 




The jet in NGC 4258 is observed at milliarcsecond 
scales in 22 GHz continuum with the VLBA ([Her 


rnstein et al. (1997a)) and at arcsecond and larger 


Evidently 10% variations in extinction can lead to 
factoy=af-two errors in the intrinsic luminosity. 


lima 


have reported X-ray 


servations of NGC 4258 with ASCA. They fit a four 
component model to the spectrum, including an ob¬ 
scured power law, an iron line, a Raymond-Smith 
plasma, and thermal bremsstrahlung. They find an 
integrated luminosity for the power law component 
in the 2 to lOkeV band of (4.4 ± 1.1) x 10 4 °ergs -1 
(at 7.3Mpc), and an absorbing column Njj = (1.5 ± 
0.2) x 10 23 enr -2 Zq/Z. Assuming a cosmic abundance 
of helium, this implies £ = 0.36 genr -2 (Zq/Z). The 
X-ray column is a factor of 4 larger than the col¬ 
umn inferred from infrared extinction which, given 
the uncertainties, may be considered reasonably good 
lent Tnripprl flip two minimis need not agree 


agreepi 

if thd IR emitting material is also producing the X- 

ray absorption. The integrated luminosity implies 
vL v =j= 2 7 x 10 40 ergs -1 at 4 5kcV, taking the mean 


of v By across the 2 to lOkeV band . 


Makishima et al.’s data are now publicly available 
through the ASCA data archives. We have tried 
fitting only an absorbed power law to the 2.5 to 
lOkeV data. While the total luminosity is robust, 
being mostly determined by the flux near lOkeV, our 
absorption column is lower than the Makishima et 
al. v alue _ This suggests that the derived absorption 


va. 

mL 


column is sensitive to the presence of the thermal 


bremsstrahlung component in the Makishima et al. 
model, and may be somewhat less reliable than the 
total luminosity. 

Herrnstein et at. (1998) report a 3cr upper limit on 
the 22 GHz flux from the neighborhood of the com- 
pact object of 0.22 mJy, or vL u = 3.1 x 10 35 ergs -1 . 

Rieke &: Lebosky (1978) report a detection in a 
5" aperture at lO^rn of 100 mJy, or vL v = 1.9 x 
10 41 ergs -1 . A consistent result at similar aperture 
size is reported by pizdziel et al. (1985). 


Because of 

the large aperture these are upper limits. 

Finally, the nucleus has been detected in polarized 
optical light by Wilkes et al. (19951 The polar- 
ized continuum has /,, ~ i/ _J - 1±u , with an amplitude 
zzLj, = 1.1 x 10 39 ergs -1 . The polarized spectral lines 
have FWHM of order 2500 km sec -1 , a velocity char¬ 
acteristic of radii just within the radius of the mas- 
ing disk. These results offer only weak constraints on 
M because the composition and distribution of the 
scattering medium are poorly constrained. Following 
Wilkes et al. we will regard the polarized emission as 


providing an upper limit on of 4 x 10 43 


3.2. Jet 


ergs 


-i 


scales in the radio, optical, and X-ray. 

The VLBA jet was discovered by [Herrnstein et 


al. (1997a)|, while higher signal-to-noise data has re 


cently been presented by Herrnstein et al. (1998) 


The northern jet has flux density ~ 3mJy, is vari¬ 
able on a timescale of weeks, and is unresolved. Its 
mean location is about 0.014 pc (= 17 light-days 
= 7500GM/c 2 ) north of the best-fit location for 
the compact object, and its position varies signifi¬ 
cantly with time. The southern jet has flux den¬ 
sity ~ 0.5 mJy, has not varied significantly in flux 
or position, and is also unresolved. It is located 
about 0.035pc (= 43 light-days = 19,000 GM/c 2 ) 
south of the best-fit location for the compact ob¬ 
ject. The difference in brightness between northern 


and southern jet is attributed by Herrnstein et al 


(1997a) to free-free absorption in the rnasing disk. 


The emitting blobs are approximately due south and 


north of the compact central object ( Herrnstein et al 
(1998)| ). This is consistent with the jet being normal 


to the outer, rnasing disk and suggests that the disk 
is not strongly warped, at least in position angle, be¬ 
tween the rnasing disk and the radius where the jet 
is launched. 

Radio and X-ray observations show a curious, 
twisted jet outflow both to the north and the south 
extending out to ~ 5 kpc (summarized in Cecil et 


al. (1992) , peril et al. (1995a)| , peril et al. (19951)) ) 
The jet X-ray emission has a power ~ 2 x 10 4U erg s 1 
and can be satisfactorily fit by a thermal spectrum 
with temperature Tx ~ 3 x 10 6 K. Associated with 
this jet are forbidden emission lines which appear to 
come from locally photoionized shocked plasma. The 
line widths suggest velocities ~ 500 km s -1 . 

4. INNER DISK MODELS 

An observed flux F 0 f, s suggests a minimum accre¬ 
tion rate M m in = 47r d 2 F 0 b s / (ec 2 ), where e = 0.1e_i is 
the accretion efficiency. This is not a rigorous limit 
since the flow can emit anisotropically and a spin¬ 
ning black hole can provide an additional reservoir 
of energy. Setting aside these concerns for the mo¬ 
ment, it is reasonable to take e_i ~ 1 (recall that 
the specific binding energy on the last stable circular 
orbit varies from 0.038c 2 for retrograde orbits around 
a maximally rotating black hole to 0.42c 2 for pro¬ 
grade orbits around a maximally rotating hole; this 
binding energy is approximately related to the effi¬ 
ciency e of a thin disk accretion flow). Taking the 
observed flux densities in the infrared, integrating 
over the band, and adding this to the 2 — lOkeV 
flux (while this is corrected for absorption, the to- 
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tal flux is not sensitive to the absorption), we find 
M m in = 10 - 5 e_i -1 Mq yr -1 . This assumes the rest 
of the spectrum is black (which is not correct given 
the detection of scattered light in the optical by 
) and the infrared extinction is 
zero. A less radically skeptical approach is to set the 
luminosity equal to the extinction-corrected infrared 
vLv. Then M m i n = 7 x 10 - 5 (2.1) s e_i -1 Mq yr -1 . 

The rest of the spectrum is not black. By analogy 
with other AGN, NGC 4258 should have a rather flat 
spectrum in uL v . One can then ask what is the mini¬ 
mum luminosity in a power law spectral energy distri- 
butio ET~that extends from the optical to 100 keV and 

passes through the X-ray data point; we find the min¬ 
imal power law slope is vL v ~ zA- 27 , and the corre¬ 
sponding luminosity is 2.3 x 10 41 ergs -1 . Adding this 
to the minimal accretion rate from the infrared point 
gives Mmin = (7(2.l) s + 4) x 10 - 5 e_i -1 M e yr -1 , al¬ 
though some of the optical, UV, and soft X-ray flux 
may be absorbed by a disk and reradiated in the in¬ 
frared. 

So far we have assumed that the accretion flow 
radiates isotropically, but a thin disk, for example, 
has f v ~ cos(z), if we ignore the effects of limb- 
brightening. Since the masing disk inclination is 82°, 
and the most natural assumption is that the inner 
disk has the same inclination (although there is no 
direct evidence for this), there will also be a large ad¬ 
ditional upward correction in the mass accretion rate, 
since cos(82°) ~ 0.14. A specific model is required for 
a quantitative estimate of the inclination correction. 

4.1. Thin Disk Model 

Suppose that the central engine is just a thin, flat 
disk that radiates locally like a black body. On the 
zA / 3 portion of the classical thin disk spectrum, 

f v = 17.8M^4 3 cos(i)A -1 / 3 mJy, ( 8 ) 

where A is in microns and M _4 is the accretion rate 
in units of 10 -4 Mq yr -1 . Inverting this equation, the 
required mass accretion rate is 

M ~ 1.3 x 10 - 6 / 3 / 2 A 1/2 cos(z ) -3/2 M 0 yr -1 , (9) 

where f u is in mJy. Assuming f u ~ zA/ 3 and 
i = 82°, the infrared data of Chary &; Becklin require 
M = 0.011 Mgyr -1 . Notice that the inclination cor¬ 
rection is quite large because M depends nonlinearly 
on cos(i). This also assumes the disk extends down 
to r = 0 ; if the disk is truncated properly (say at 
r = 6 GM/c 2 ) then M must be evaluated numeri¬ 
cally: M = 0.013 Mq yr -1 . Thin disk spectra for 
i = 0°,45°, and 82° are shown in Figure 1. 


Wilkes et al. (1995) 


There are several effects that might tend to change 
the accretion rate. First, the inner disk might be in¬ 
clined differently from the masing disk. If i = 45°, 
then the accretion rate can be pushed as low as 
0.0012 Mq yr -4 ; the order of magnitude change in M 
again follows from the nonlinear dependence of M 
on cos(i) in eq.(j|). It is also possible, however, that 
the inclination of the inner disk is larger than 82°. 
Second, the obscuring dust grains might well have a 
different size distribution than in the local interstellar 
medium. The sense of the effect of grain size evolu¬ 
tion, however, at least in the grain models of Laor & 
Draine (1993)|, is to decrease Ah/Ak and so increase 


the required mass accretion rate. 

Third, some of the near-infrared emission may be 
reprocessed light from the central engine, either in 
a disk or in a wind (e.g. [Konigl &; Kartje (1994) ). 
Crudely speaking, the light must be reprocessed at 
risoo ~ 2.8 x ltfl}£GM/c 2 corresponding to a 


f r*j 


characteristic temperature of 1500A" for near IR emis¬ 
sion (dust grains much smaller than 1 /zrn can produce 
H and K band emission from a somewhat larger ra¬ 
dius because they radiate inefficiently in the infrared). 
When is reprocessing important in a disk? Suppose 
/n(r, cj)) is the “obliquity” of the disk, which is to say, 
the angle between the normal to the disk surface and 
the radius vector, / is the fraction of intercepted ra¬ 
diation that is thermalized in the disk (so 1 — / is 
an effective albedo). Then reprocessing is dominant 
if r A (GM/c 2 )(e/cos^z) -1 . For a thin, flared disk 
cos /z ~ H/R. One can then show that, in a standard 
a disk model over a large range of M, reprocessing is 
not dominant at r ~ risoo- 

If the disk is warped then cos /z can be larger and 
reprocessing can dominate internal heating. It is not 
possible to survey the warp parameter space, since 
a warp is described by two free functions: the incli¬ 
nation i(r) and position angle p(r). A broad range 
of behavior is possible. For example, the warp could 
be made sufficiently large so as to cover the central 
engine (raising the inferred accretion rate); it could 
also be arranged so that the inner disk is nearly face 
on (lowering the inferred accretion rate). 

As an illustration, however, we have constructed a 
“naive” warped disk model with only an inclination 
warp (p(r) = const.) in which the reprocessed light 
is reradiated as a black body. The disk extends from 
r = 6GM/c 2 to 0.26 pc, as does the flat disk model, 
and has constant M. The effective temperature T e of 
the disk is given by T e 4 = Tf nt + fobiTexti where T int is 
the effective temperature in the absence of external 
heating, and 


°T* xt = 


eMc 2 
47 T r 2 


( 10 ) 
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and e = 0.1. Here f 0 u is the obliquity factor. If 
the disk is illuminated on the side facing us, then 
fobi — sin((j))\di/dlnr\, where 4> is the disk orbital 


An alternative to the thin disk plus corona model 
for the inner accretion flow is the two-temperature 
advection-dominated accretion flow model (ADAF; 


phase. If the disk is illuminated on the side facing 

Ichimaru (1977), Narayan & Yi (1994), |Narayan & 

away from us, as it is for half its orbit, then the re- 

Yi (1995b), 

Abramowicz (1995)|; see|Narayan (1997), 

processing only adds a flux <jT ^ xt /r to our side of the 

|Narayan et al. (1998) for reviews). ADAFs are opti- 

disk, and so, assuming that r 3> (T ext /T int ) 4 , we set 
fobl = 0. 

cally thin and geometrically thick, with H/R ~ 1. 
They have been used to successfully fit the spec- 


Figure 2 shows the resulting thin, warped disk spec¬ 
trum. The warp extends from r = 100 GM/c 2 to 
r = 0.26pc with di/dlnr = 0.05. The accretion rate 
is 8 x 10” 4 Mq yr^ 1 . Three effects combine here to re¬ 
duce the required accretion rate: (1) reprocessing en¬ 
hances the 2/rin flux; (2) reprocessing make the spec¬ 
trum shallower at 2/im thus reducing the reddening 
correction; (3) the inclination warp changes the incli¬ 
nation of the inner disk from 82° to 61°, thus making 
the reprocessing portion of the disk more face-on. 

Our warped disk model is naive in the sense that 
the reprocessed emission does not emerge as a black- 
body. In particular, dust grains in the disk atmo¬ 
sphere can cause half of the flux intercepted by the 
disk to reemerge as a dilute blackbody with T > T e . 
This effect has been considered in the context of 
quasar disks by Sanders et al. (1989)| and Phinney 


(19890; the same effect has been considered in the 
context of circumstellar disks by |Chiang fc Goldres) 
ich (1 |997) . In Phinney’s quasar models the resulting 
spectrum has a notch at ~ 1/rrn related to the onset 
of dust sublimation at T ~ 1800 K. Since the spec¬ 


tral energy distribution of accretion flows onto both 
stellar mass and supermassive black holes. In their 
simplest form, ADAFs approximate the flow as non- 
radiative, with nested spherical surfaces moving ra¬ 
dially inwards with an angular velocity given by the 
value on the equatorial plane. Radial pressure gradi¬ 
ents combine with centrifugal force to oppose grav¬ 
ity and radial motion is fast, though subsonic until 
it approaches the event horizon. There is hypothe¬ 
sized to be a large viscous stress (~ aP, where P 
is the pressure) that transports angular momentum 
radially outward. 

The ADAF approach has some dynamical short¬ 
comings. Firstly, it does not take account of an¬ 
gular velocity gradients on spherical surfaces. This 
appears to be a reasonable approximation far from 
the black hole in the self-similar region of the flow 
(cf. Narayan fe Yi (1995a)|), but the approxima¬ 


tion may break down close to the hole, especially 
close to the symmetry axis. This is relevant for un¬ 
derstanding whether the flow forms a funnel and a 
jet. The indications are that there is no funnel when 


trum might then be declining at H and K bands, this 

« 0.01 ( 

Narayan, Kato & Honma (1997), Gam- 

could reduce the reddening correction and thus also 

mie &; Popham (1998), 

Popham &; Garnmie (1998), 

reduce the required accretion rate. 

though this conclusion 

is not certain). Secondly, 


Reprocessing is ultimately limited, however, by the 
luminosity of the central accretion flow. Figure 1 also 
shows a thin disk spectrum for a face-on disk trun¬ 
cated at r = 6 GM/c 2 , with M = 7 x 10~ 6 MQyr~ 1 , 
corresponding to the model of Neufeld &; Maloney! 


(1995[j] with a = 0.1. It would be difficult to produce 
the observed infrared flux from such a low accretion 
rate ( |lisk even with a highly optimized arrangement 
of reprocessing material. 

A reprocessing model for the H and K band flux 
can be tested via long term monitoring. If the light 
is reprocessed, then IR variability should be highly 
correlated with X-ray variability, with delays on or¬ 
der the light travel time nsoo/c — 5.5 x 10 5 L 1 J 2 ~ s. 

Finally, while we have not explicitly fit the X- 
ray emission using a thin disk plus corona model 
(e.g. Haardt Sz Maraschi (199ll) , it is generally ac¬ 
knowledged that it is possible to do so, since there 
are a number of free parameters associated with the 
corona. 

4.2. Advection-Dominated Flow Model 


the solutions assume vertical hydrostatic equilibrium, 
which is suspect in the supersonic region close to 
the hole. This approximation should not have di¬ 
rect impact on the emergent radiation, though. Fi¬ 
nally, in these solutions, the gas has positive Bernoulli 
constant and always has enough e n ergy to escape 
the hole (see Narayan &; Yi (1994) , Narayan &; Yi 
(1995a)| , Narayan, Kato &; Honma (1997)| )) This 


arises because of the inevitable, outward radial trans¬ 
port of energy by the viscous stress, which is balanced 
exactly by inward advective transport of energy. This 
balance must be maintained everywhere. If it is not, 
then it is possible that energy may be carried away 
by a wind (see Blandford (1998)| , and in preparation). 

Counterbalancing these shortcomings is the fact 
that the ADAF is the only model of hot accretion 
flows that attempts to calculate the dynamics, ther¬ 
mal structure and radiative spectrum of the gas self- 
consistently. Other approaches, e.g. coronal models, 
leave most details of the dynamics unspecified and 
allow far greater freedom in choosing the properties 
of the radiating plasma. 
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Lasota et al. (1996) proposed a model for NGC 
4258 which comprises a thin disk extending outward 
from the transition radius ?’t r and an ADAF extend¬ 
ing inward to the black hole. The mass accretion 
rate in the model is essentially fixed by the X-ray 
data point to be of order 0.014a_i Lasota 

et al. take a_i = 1. One free parameter remains: 
rtr- At the time Lasota et al. (1996) was written the 
data did not constrain r tr . Two results have since ap- 


peared that do: the 22 GHz upper limit of Herrnstein 
et al. (1998), and the IR detections of Chary fe Beck 


lin (1 197)| . There has also been considerable improve¬ 
ment in the calculation of the ADAF model spectral 
energy distribution. Relativistic effects (|Gammie &| 
Pophmi (1998), and references therein) and advec- 
tive heating of the electrons ( Nakamura et al. (1997) ) 
are now included; details of the computation are de¬ 
scribed in Narayan et al. (1998) , which applies the 
ADAF model to the Galactic Center source Sgr A*. 

Figure 3 shows the spectral energy distribution 
of NGC 4258 corresponding to three ADAF mod¬ 
els, with = (10,10 1 ' 5 ,10 2 ) X 2GM/c 2 and M = 
(0.012,0.008, 0.006) yr -1 , which we call models 

1,1.5, and 2. The disk outside rt r is thin, flat, and ra¬ 
diates like a blackbody. All models use the standard 
values for other parameters: a = 0.3 , (3 = 0.5, and 
5 = ICE 3 (f3 = Pgas/Ptot , 5 = fraction of dissipation 
going into electrons). Model 1 has the most (least) 
flux in the optical (radio). Evidently model 1 fits the 
IR data best and does not violate the 22 GHz upper 
limit, although the model 1.5 also gives an acceptable 
fit to the IR data and may have better agreement 
with the spectral index of the scattered optical light. 
Model 2 is close to violating the 22 GHz upper limit. 

The ADAF model predicts the flux from the cen¬ 
tral accretion flow at several wavelengths that may 
be observed in the near future. At 10/rm the ob¬ 
served flux should not be below the thin disk value 
(see equation J|]) of (200,150,130) cos (i) mJy for 
models (1,1.5, 2); it can be larger because of rera¬ 
diation by dust. At 22 GHz the flux should be 
(0.017,0.082,0.23) mJy with s = din ///din u ss 
1.5 for all three models. Higher frequency in¬ 
terferometry may be possible; at 100 GHz, f u = 
(0.13,0.65,1.5) mJy, with s = (1.8,1.4,1.3), and at 
690 GHz,/,, = (3.4,10., 18.) mJy, s = (1.6,1.3,1.0). 
At 5500A the expected flux is (38,8.4,0.59) mJy with 


s = (—0.39,—2.0,—2.2); recall that Wilkes et al 
(1995[| find s = —1.1. The model is adjusted so 
that the flux at lOkeV fits the observations; at 
lOkeV the predicted photon index is (1.9,1.9, 2.0), 
which may be compared with the observed index of 
1.78 ± 0.29. Finally, the expected flux at 30keV is 
(3.1,3.4, 3.4) x 10~ 6 phcm~ 2 s^ 1 keV” 1 , and the pho¬ 
ton index for all three models is 2.0. 


4.3. Ion Torus Model 

A related, though distinct, alternative to the ADAF 


is the ion torus ( Reesetal. (1982) ). As in the ADAF, 
it is assumed that the plasma is two-temperature with 
most of the pressu re being supplied by the ions (see 
Shapiro et al 1976 for the first application of the two- 
temperature idea to accretion) and that the radia¬ 
tive losses are small. The viscosity is assumed to be 
small enough (low a), however, that poloidal motion 
can be ignored in the equation of hydrostatic equilib¬ 
rium. In its simplest manifestation, the specific an¬ 
gular momentum is constant and the surface of the 
torus is a surface of constant binding energy. As the 
radiative losses are small, the binding energy is quite 
low (relative to c 2 per unit mass) and plasma flows 
onto the black hole through an equatorial cusp with 
this energy and consistent with having lost little of 
its binding energy. This approach can accommodate 
the peculiarly relativistic features of the Kerr metric 
( Abramowicz, Jaroszyriski fe Sikora (1978)| ). 

Following previous work on thick disks, the ion 
torus model simply specifies the angular momentum 
profile rather than solving an angular momentum 
equation, and this is a serious weakness. Also, the 
model has no ready prescription for describing the 
gas inflow and, consequently, the thermal structure. 
In addition, tori of this sort are subject to hydro- 
dynamical instability ( |Papaloizou fe Pringle (1984) ), 
which calls into question the assumption of slow an¬ 
gular momentum transport. 

Nonetheless, ion tori do have one feature that is 
important in the present context: they naturally cre¬ 
ate a pair of funnels in which magnetic flux can be 
concentrated so as to extract energy from the cen¬ 
tral spinning black hole and power the jets. Accord¬ 
ing to the work of Narayan, Kato &; Honma (1997) 


and Pammie fc Popliam (1998) , funnels are likely 
to be most prominent in low-a flows (ion tori), and 
are much less pronounced in high-a flows (standard 
ADAFs). 


4.4. General Constraints 

The failure to detect radio emission from the loca¬ 
tion of the black hole has some more general implica¬ 
tions. Although we have demonstrated that ADAF 
models can be made to be consistent with this upper 
limit, alternative models that admit a non-thermal 
distribution of relativistic electrons to be accelerated 
directly are quite tightly constrained. 

Let us suppose that there is a radio photo¬ 
sphere where the optical depth to synchrotron self¬ 
absorption becomes unity. Let it be located at a ra¬ 
dial distance r p = 10 14 r p i 4 cm from the black hole. 
The 22 GHZ upper bound on the flux implies that 
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the radio brightness temperature of this photosphere 
must satisfy 

T < 2 x lO 11 ^ 2 K (11) 

flHerrnstein et al. (1997a)| ). The corresponding en¬ 
ergy of the electrons emitting at 22 GHz is ~ 3kT ~ 


E p ~ 


70 B 


- 1/2 


MeV. These electrons radiate at 


7 2 B MHz (cf. Rybicki fc Lightman (19791) , and 
so we conclude that the magnetic field strength at 
the photosphere satisfies B p = B(jr p ) ^ r 4 14 G, (as¬ 
suming that it is significantly less than the cyclotron 
field, ~ 7000 G). In other words the magnetic pres¬ 
sure must satisfy P mp ^ 0.04r 3 14 , measured in dyne 
cm -2 . We use the theory of synchrotron radiation to 
solve for the radius of the 22 GHz photosphere: 


r P i4 ~ 10 f e P; 


2 

mp’ 


( 12 ) 


where f e is the ratio of the emitting electrons to 
the magnetic pressure. Now, in order not to ex¬ 
ceed the radio flux upper limit, we require that 
r p i 4 <; 1.3/g' 07 ~ 20(GAf/c 2 )/g' 07 and, equivalently, 
P mp <; 1.1 /(T 0 03 dyne cm” 2 . As emphasized below, 
however, typical pressures in accretion flows are gen¬ 
erally expected to exceed ~ 1 dyne cm” 2 by a large 
factor, so this is a fairly stringent constraint. 

ADAF models satisfy this constraint because they 
ignore extraneous relativistic particle acceleration so 
that the electron distribution function is thermal and 
exponentially declining at energies below E p . (Ma- 


liade fan &; Quataert (1997) have shown that the ther- 
malization time scale for the electrons in an ADAF is 
shorter than the flow time scale, so that the assump¬ 
tion of a thermal particle distribution is consistent.) 
What we see from the analysis in this section is that 
this assumption is necessary, - 

Bu lk corona models generally do postulate rolativis 

tic particle acceleration. The accelerated particles 
might be expected to cause observable radio emission, 
and so the limit on the 22 GHz brightness tempera¬ 
ture is a serious constraint. 

There are some processes which can suppress radio 
synchrotron emission. The Razin effect is important 
at densities above ~ 10 9 B p cm” 3 . Secondly, relativis¬ 
tic Coulomb interactions and synchrotron loss may 
cool freshly accelerated electrons to keep their steady 
pressure at a low level. This happens on time scales 
competitive with the inflow timescale. Thirdly, ex¬ 
trinsic free-free absorption may extinguish the radio 
source. Two difficulties with this explanation are that 
the southern jet can be seen through the accreted gas 
at a projected perpendicular distance from the black 
hole of ~ 6000GM/c 2 and that there is no observed 
soft X-ray absorption. 

5. JET 


It is clearly of interest to relate the observed VLBI 
jets to the gas flow close to the black hole. If we con¬ 
sider the northern VLBI component, we find that, at 
a projected distance of rj_ ~ 4 x 10 16 cm, the 22 GHz 
flux density is ~ 3 mJy. If we suppose that the unre¬ 
solved source is the base of a jet with an opening angle 
of ~ 10° (cf. Herrnstein et al. (1997a) (propagating 


nearly perpendicular to the line of sight with a bulk 
Lorentz factor T, then the radio brightness temper¬ 
ature is only ~ 2 x 10 9 K. Even allowing for rela¬ 
tivistic beaming this is too low for synchrotron self¬ 
absorption to be important and it seems reasonable to 
assume that both radio components are optically thin 
at this frequency. Let us now assume that the rela¬ 
tivistic energy density in the jet is not much smaller 
than the magnetic energy density. A straightforward 
calculation relates the minimum power carried by the 
relativistic electrons (and, quite possibly, positrons) 
in the jet to T and to the ratio of the electron energy 
density to the magnetic energy density /3 e (assumed 
> 1 ): 


L jet ~ 7 x 10 42 /3 9 ' 4 


3.4 


erg s 


-i 


(13) 


(see, e.g., Blandford (1990)). The strong dependence 
of the minimum jet power upon the Lorentz factor 
is due to a combination of the bulk energy and the 
diminution of the observed flux by relativistic beam¬ 
ing when the observer is located roughly perpendic¬ 
ular to the jet direction. We see that if jets have ul- 
trarelativistic flow speeds, then they probably carry 
off significantly more power than is radiated directly. 

The power in the jet may be extracted directly 
from the spin of the black hole (Blandford & Znajek 


(1977)| ). If instead, the power is derived from the ac¬ 


cretion flow, then the jet implies a certain minimum 
mass accretion rate 


M ~ 1.2 x 10” 3 e_i” 1 /?g' 4 


3.4 


M© yr” 1 . (14) 


Notice that this is similar to the estimates of M ob¬ 
tained earlier by fitting the observed infrared spec¬ 
trum. 

The pressure in the jet ~ 0.07/3g' 4 (r/10) L4 (r/4 x 
10 16 cm)” 2 dyne cm” 2 must be confined transversely. 
If we evaluate this at the transition radius of 
the ADAF solution, then the jet pressure is ~ 
700(r/10) 1 ' 4 /3g' 4 dyne cm” 2 . For comparison, the 
pressure in the ADAF flow is 


P ~ 8000 


M 


10” 2 M©yr 


-l 


a- 


rT 


Hr 


V60 GM/c 2 


-5/2 


dyne cm 
(15) 



















We i.ote that, although the formation of funnels _ &; Maloney (1995) from consideration of the masing 


appears to be problematic in ADAFs with a ^ 
0.01 ( Narayan, Kato fc Honma (1997) , pammie~~& 
Poph rm (1998)| ), the pressure required to confine a 
jet able to account for the observed radio sources is 
probably present. An ion torus, by contrast, natu¬ 
rally forms a torus and can exert more pressure than 
a high-a ADAF as long as it is unable to cool. If 
the accretion rate is 10 -4 M 0 yr~ 4 , however, it 
does not seem possible to supply enough pressure to 
confine a relativistic jet. 


disk. 

The minimum power in the jet in NGC 4258 also 
constrains the mass accretion rate, provided the jet 
power is derived from the accretion flow and not from 
the spin of the black hole. The resulting estimate for 
M (cf. equation |l4]]) is similar to that obtained by 
fitting the infrared spectrum. 

The limit on the 22 GHz radio flux in NGC 4258 
constrains the transition radius in any ADAF to less 


than 100 Schwarzschild radii, as emphasized by Her 


rnstein et al. (1998). As we show in this paper, 


6. DISCUSSION 

We have argued that the outer, masing disk in 
NGC 4258 is unlikely to be in a steady state because 
of its long viscous timescale. Thus it is questionable 
whether conditions in the masing disk can be used 
to estimate the mass accretion rate onto the central 
black hole. 

We have fit several models to the spectral energy 
distribution of the nucleus in order to better con¬ 
strain the accretion rate. The ADAF model is able 
to fit all the data as long as M ~ 10~ 2 MQyr _1 , 
and the transition radius (where the accretion flow 
makes the transition from a thin disk to an ADAF) 
is ~ (10 — 100 )GM/c 2 . The ADAF accretion rate 
is basically fixed by the X-ray luminosity, but it is 
nicely consistent with the observed infrared flux as 
long as the infrared is generated by a thin, flat disk 
with inclination close to that of the masing disk. 

We have also fit pure thin disk models to the in¬ 
frared data. (The X-ray flux would presumably be 
generated by a hot corona, but we have not con¬ 
structed an explicit corona model, since coronae are 
not as tightly constrained as ADAFs.) The inferred 
accretion rate is strongly inclination dependent. For 
i = 82°, M = 0.013 Mq yr^ 1 . A lower accretion 
rate of 0.0012 Mq yr _1 is obtained by giving the disk 
an inclination of 45°, although such a low inclina¬ 
tion seems unlikely given the orientation of the large 
scale jet. The accretion rate can be further lowered 
by giving the disk a warp that reprocesses radiation 
from the central source into the infrared. Since the 
warp parameter space is large, we have calculated 
the spectrum of a representative model disk with a 
modest warp, and the implied accretion rate is 8 x 
10~ 4 MQyr _1 . The accretion rate might be pushed 
lower by a somewhat artificial positioning of repro¬ 
cessing gas at the right radius; this could push the 
accretion rate as low as 7 x 10 _5 (2.1) s e_i _1 Mq yr _1 , 
where s is the slope of the intrinsic spectrum. For 
s = —1, e_i = 1, and a_i = 1 (larger values seem 
unlikely except possibly in the event of gravitational 
instability), this minimal accretion rate is still larger 
than the 7 x 10 -6 a_i Mq yr _1 obtained by [Neufcld| 


the flux limit also constrains corona models. If the 
observed X-ray emission in NGC 4258 is from a 
corona, then the corona must be primarily thermal 
rather than non-thermal, since otherwise it would 
produce more radio flux than observed. This limits 
the kinds of relativistic acceleration one can invoke 
in the corona. 

New data at several different frequencies would give 
interesting information on the nature of the central 
accretion flow in NGC 4258. The ADAF model pre¬ 
dicts a complete spectrum whose amplitude and slope 
can be compared with new data in the radio, mil¬ 
limeter, infrared, and X-ray; the predicted values are 
given in §4. 

It would be particularly interesting to know if the 
near infrared flux is reprocessed or not. This can 
be tested by looking for variation on a timescale 
t r = risoo/c — 5.5 x 10 5 L42 2 s, which is about a 
week. Since the viscous timescale- indeed, even the 
dynamical timescale- is much longer than this in the 
disk at risoo, we would not expect significant varia¬ 
tion in this component of the flux unless the infrared 
flux is reprocessed, or the infrared flux is produced 
nonthermally. If a significant fraction of the H and 
K band flux is reprocessed, this would rule out our 
simplest models and we would be forced to invoke a 
warped disk, wind, or other arrangement of repro¬ 
cessing gas. Higher angular resolution mid-infrared 
data would also be interesting, since a flat disk pre¬ 
dicts a steeper spectrum in the mid-infrared than a 
warped, reprocessing disk. 

Finally VLBI observations at lower frequencies 
than 22 GHz may tighten the general constraints on 
relativistic particle acceleration close to the black hole 
and exclude free-free absorption as a factor. They 
may also be able to determine the onset of self¬ 
absorption in the jet and thence provide a stronger 
constraint on the pressure and jet power by allowing 
us to solve for the magnetic field strength. Further 
monitoring of the variability in both the northern and 
the southern component may also help determine the 
jet speed and elucidate the relationship between the 
jets and gas flow within ~ 100 gravitational radii of 











9 


the central black hole. 
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Fig. 1. The spectral energy distribution of the nucleus of NGC 4258 . The data points are described in the text (see §2) 
The models are flat, thin disks with inclination i = (45, 82, 0)° and M = (0.0012, 0.013, 7 x 10~ 6 ) M 0 yr _1 . 
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Fig. 2.— Same as Figure 1, but for a warped, reprocessing disk wit 
= 82° and M = 0.013 M 0 yr _1 . 



M = 0.0008 Mq yr 1 as well as the flat disk model with 
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Fig. 3.— The models are ADAFs fit with rtr = (10,10 1 ' 5 ,10 2 )2GAf/c 2 and M = (0.012, 0.008, 0.006) M® yr As the transition 
radius moves in the optical/IR flux increases. 



